Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures.

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

Structural and mechanistic basis of the central energy-converting methyltransferase complex of methanogenesis.

Methanogenic archaea inhabiting anaerobic environments play a crucial role in the global biogeochemical material cycle. The most universal electrogenic reaction of their methane-producing energy metabolism is catalyzed by N    5-methyl-tetrahydromethanopterin: coenzyme M methyltransferase (MtrABCDEFGH), which couples the vectorial Na+ transport with a methyl transfer between the one-carbon carriers tetrahydromethanopterin and coenzyme M via a vitamin B12 derivative (cobamide) as prosthetic group. We present the 2.08 Šcryo-EM structure of Mtr(ABCDEFG)3 composed of the central Mtr(ABFG)3 stalk symmetrically flanked by three membrane-spanning MtrCDE globes. Tetraether glycolipids visible in the map fill gaps inside the multisubunit complex. Putative coenzyme M and Na+ were identified inside or in a side-pocket of a cytoplasmic cavity formed within MtrCDE. Its bottom marks the gate of the transmembrane pore occluded in the cryo-EM map. By integrating Alphafold2 information, functionally competent MtrA-MtrH and MtrA-MtrCDE subcomplexes could be modeled and thus the methyl-tetrahydromethanopterin demethylation and coenzyme M methylation half-reactions structurally described. Methyl-transfer-driven Na+ transport is proposed to be based on a strong and weak complex between MtrCDE and MtrA carrying vitamin B12, the latter being placed at the entrance of the cytoplasmic MtrCDE cavity. Hypothetically, strongly attached methyl-cob(III)amide (His-on) carrying MtrA induces an inward-facing conformation, Na+ flux into the membrane protein center and finally coenzyme M methylation while the generated loosely attached (or detached) MtrA carrying cob(I)amide (His-off) induces an outward-facing conformation and an extracellular Na+ outflux. Methyl-cob(III)amide (His-on) is regenerated in the distant active site of the methyl-tetrahydromethanopterin binding MtrH implicating a large-scale shuttling movement of the vitamin B12-carrying domain.

Using cryo-EM to understand the assembly pathway of respiratory complex I.

Complex I (proton-pumping NADH:ubiquinone oxidoreductase) is the first component of the mitochondrial respiratory chain. In recent years, high-resolution cryo-EM studies of complex I from various species have greatly enhanced the understanding of the structure and function of this important membrane-protein complex. Less well studied is the structural basis of complex I biogenesis. The assembly of this complex of more than 40 subunits, encoded by nuclear or mitochondrial DNA, is an intricate process that requires at least 20 different assembly factors in humans. These are proteins that are transiently associated with building blocks of the complex and are involved in the assembly process, but are not part of mature complex I. Although the assembly pathways have been studied extensively, there is limited information on the structure and molecular function of the assembly factors. Here, the insights that have been gained into the assembly process using cryo-EM are reviewed.

Conformational coupling of redox-driven Na+-translocation in Vibrio cholerae NADH:quinone oxidoreductase.

In the respiratory chain, NADH oxidation is coupled to ion translocation across the membrane to build up an electrochemical gradient. In the human pathogen Vibrio cholerae, the sodium-pumping NADH:quinone oxidoreductase (Na+-NQR) generates a sodium gradient by a so far unknown mechanism. Here we show that ion pumping in Na+-NQR is driven by large conformational changes coupling electron transfer to ion translocation. We have determined a series of cryo-EM and X-ray structures of the Na+-NQR that represent snapshots of the catalytic cycle. The six subunits NqrA, B, C, D, E, and F of Na+-NQR harbor a unique set of cofactors that shuttle the electrons from NADH twice across the membrane to quinone. The redox state of a unique intramembranous [2Fe-2S] cluster orchestrates the movements of subunit NqrC, which acts as an electron transfer switch. We propose that this switching movement controls the release of Na+ from a binding site localized in subunit NqrB.

Cyclic di-AMP traps proton-coupled K+ transporters of the KUP family in an inward-occluded conformation.

Cyclic di-AMP is the only known essential second messenger in bacteria and archaea, regulating different proteins indispensable for numerous physiological processes. In particular, it controls various potassium and osmolyte transporters involved in osmoregulation. In Bacillus subtilis, the K+/H+ symporter KimA of the KUP family is inactivated by c-di-AMP. KimA sustains survival at potassium limitation at low external pH by mediating potassium ion uptake. However, at elevated intracellular K+ concentrations, further K+ accumulation would be toxic. In this study, we reveal the molecular basis of how c-di-AMP binding inhibits KimA. We report cryo-EM structures of KimA with bound c-di-AMP in detergent solution and reconstituted in amphipols. By combining structural data with functional assays and molecular dynamics simulations we reveal how c-di-AMP modulates transport. We show that an intracellular loop in the transmembrane domain interacts with c-di-AMP bound to the adjacent cytosolic domain. This reduces the mobility of transmembrane helices at the cytosolic side of the K+ binding site and therefore traps KimA in an inward-occluded conformation.

Insights into complex I assembly: Function of NDUFAF1 and a link with cardiolipin remodeling.

Respiratory complex I is a ~1-MDa proton pump in mitochondria. Its structure has been revealed in great detail, but the structural basis of its assembly, in humans involving at least 15 assembly factors, is essentially unknown. We determined cryo-electron microscopy structures of assembly intermediates associated with assembly factor NDUFAF1 in a yeast model system. Subunits ND2 and NDUFC2 together with assembly factors NDUFAF1 and CIA84 form the nucleation point of the NDUFAF1-dependent assembly pathway. Unexpectedly, the cardiolipin remodeling enzyme tafazzin is an integral component of this core complex. In a later intermediate, all 12 subunits of the proximal proton pump module have assembled. NDUFAF1 locks the central ND3 subunit in an assembly-competent conformation, and major rearrangements of central subunits are required for complex I maturation.

The molybdenum storage protein forms and deposits distinct polynuclear tungsten oxygen aggregates.

Some N2-fixing bacteria store Mo to maintain the formation of the vital FeMo-cofactor dependent nitrogenase under Mo depleting conditions. The Mo storage protein (MoSto), developed for this purpose, has the unique capability to compactly deposit molybdate as polyoxometalate (POM) clusters in a (αß)3 hexameric cage; the same occurs with the physicochemically related tungstate. To explore the structural diversity of W-based POM clusters, MoSto loaded under different conditions with tungstate and two site-specifically modified MoSto variants were structurally characterized by X-ray crystallography or single-particle cryo-EM. The MoSto cage contains five major locations for POM clusters occupied among others by heptanuclear, Keggin ion and even Dawson-like species also found in bulk solvent under defined conditions. We found both lacunary derivatives of these archetypical POM clusters with missing WOx units at positions exposed to bulk solvent and expanded derivatives with additional WOx units next to protecting polypeptide segments or other POM clusters. The cryo-EM map, unexpectedly, reveals a POM cluster in the cage center anchored to the wall by a WOx linker. Interestingly, distinct POM cluster structures can originate from identical, highly occupied core fragments of three to seven WOx units that partly correspond to those found in MoSto loaded with molybdate. These core fragments are firmly bound to the complementary protein template in contrast to the more variable, less occupied residual parts of the visible POM clusters. Due to their higher stability, W-based POM clusters are, on average, larger and more diverse than their Mo-based counterparts.

High-resolution structure and dynamics of mitochondrial complex I-Insights into the proton pumping mechanism.

Mitochondrial NADH:ubiquinone oxidoreductase (complex I) is a 1-MDa membrane protein complex with a central role in energy metabolism. Redox-driven proton translocation by complex I contributes substantially to the proton motive force that drives ATP synthase. Several structures of complex I from bacteria and mitochondria have been determined, but its catalytic mechanism has remained controversial. We here present the cryo-EM structure of complex I from Yarrowia lipolytica at 2.1-Å resolution, which reveals the positions of more than 1600 protein-bound water molecules, of which ~100 are located in putative proton translocation pathways. Another structure of the same complex under steady-state activity conditions at 3.4-Å resolution indicates conformational transitions that we associate with proton injection into the central hydrophilic axis. By combining high-resolution structural data with site-directed mutagenesis and large-scale molecular dynamic simulations, we define details of the proton translocation pathways and offer insights into the redox-coupled proton pumping mechanism of complex I.

A hybrid approach reveals the allosteric regulation of GTP cyclohydrolase I.

Guanosine triphosphate (GTP) cyclohydrolase I (GCH1) catalyzes the conversion of GTP to dihydroneopterin triphosphate (H2NTP), the initiating step in the biosynthesis of tetrahydrobiopterin (BH4). Besides other roles, BH4 functions as cofactor in neurotransmitter biosynthesis. The BH4 biosynthetic pathway and GCH1 have been identified as promising targets to treat pain disorders in patients. The function of mammalian GCH1s is regulated by a metabolic sensing mechanism involving a regulator protein, GCH1 feedback regulatory protein (GFRP). GFRP binds to GCH1 to form inhibited or activated complexes dependent on availability of cofactor ligands, BH4 and phenylalanine, respectively. We determined high-resolution structures of human GCH1-GFRP complexes by cryoelectron microscopy (cryo-EM). Cryo-EM revealed structural flexibility of specific and relevant surface lining loops, which previously was not detected by X-ray crystallography due to crystal packing effects. Further, we studied allosteric regulation of isolated GCH1 by X-ray crystallography. Using the combined structural information, we are able to obtain a comprehensive picture of the mechanism of allosteric regulation. Local rearrangements in the allosteric pocket upon BH4 binding result in drastic changes in the quaternary structure of the enzyme, leading to a more compact, tense form of the inhibited protein, and translocate to the active site, leading to an open, more flexible structure of its surroundings. Inhibition of the enzymatic activity is not a result of hindrance of substrate binding, but rather a consequence of accelerated substrate binding kinetics as shown by saturation transfer difference NMR (STD-NMR) and site-directed mutagenesis. We propose a dissociation rate controlled mechanism of allosteric, noncompetitive inhibition.

The resolution revolution in cryoEM requires high-quality sample preparation: a rapid pipeline to a high-resolution map of yeast fatty acid synthase.

Single-particle electron cryo-microscopy (cryoEM) has undergone a 'resolution revolution' that makes it possible to characterize megadalton (MDa) complexes at atomic resolution without crystals. To fully exploit the new opportunities in molecular microscopy, new procedures for the cloning, expression and purification of macromolecular complexes need to be explored. Macromolecular assemblies are often unstable, and invasive construct design or inadequate purification conditions and sample-preparation methods can result in disassembly or denaturation. The structure of the 2.6 MDa yeast fatty acid synthase (FAS) has been studied by electron microscopy since the 1960s. Here, a new, streamlined protocol for the rapid production of purified yeast FAS for structure determination by high-resolution cryoEM is reported. Together with a companion protocol for preparing cryoEM specimens on a hydrophilized graphene layer, the new protocol yielded a 3.1 Šresolution map of yeast FAS from 15 000 automatically picked particles within a day. The high map quality enabled a complete atomic model of an intact fungal FAS to be built.

Respiratory complex I - structure, mechanism and evolution.

Respiratory complex I is an intricate multi-subunit membrane protein with a central function in aerobic energy metabolism. During the last years, structures of mitochondrial complex I and respiratory supercomplexes were determined by cryo-EM at increasing resolution. Structural and computational studies have shed light on the dynamics of proton translocation pathways, the interaction of complex I with lipids and the unusual access pathway of ubiquinone to the active site. Recent advances in understanding complex I function include characterization of specific conformational changes that are critical for proton pumping. Cryo-EM structures of the NADH dehydrogenase-like (NDH) complex of photosynthesis and a bacterial membrane bound hydrogenase (MBH) have provided a broader perspective on the complex I superfamily.

Structural basis of proton-coupled potassium transport in the KUP family.

Potassium homeostasis is vital for all organisms, but is challenging in single-celled organisms like bacteria and yeast and immobile organisms like plants that constantly need to adapt to changing external conditions. KUP transporters facilitate potassium uptake by the co-transport of protons. Here, we uncover the molecular basis for transport in this widely distributed family. We identify the potassium importer KimA from Bacillus subtilis as a member of the KUP family, demonstrate that it functions as a K+/H+ symporter and report a 3.7 Å cryo-EM structure of the KimA homodimer in an inward-occluded, trans-inhibited conformation. By introducing point mutations, we identify key residues for potassium and proton binding, which are conserved among other KUP proteins.

Analysis of the co-translational assembly of the fungal fatty acid synthase (FAS).

The yeast fatty acid synthase (FAS) is a barrel-shaped 2.6 MDa complex. Upon barrel-formation, two multidomain subunits, each more than 200 kDa large, intertwine to form a heterododecameric complex that buries 170,000 Å2 of protein surface. In spite of the rich knowledge about yeast FAS in structure and function, its assembly remained elusive until recently, when co-translational interaction of the ß-subunit with the nascent α-subunit was found to initiate assembly. Here, we characterize the co-translational assembly of yeast FAS at a molecular level. We show that the co-translationally formed interface is sensitive to subtle perturbations, so that the exchange of two amino acids located in the emerging interface can prevent assembly. On the other hand, assembly can also be initiated via the co-translational interaction of the subunits at other sites, which implies that this process is not strictly site or sequence specific. We further highlight additional steps in the biogenesis of yeast FAS, as the formation of a dimeric subunit that orchestrates complex formation and acts as platform for post-translational phosphopantetheinylation. The presented data supports the understanding of the recently discovered prevalence of eukaryotic complexes for co-translational assembly, and is valuable for further harnessing FAS in the biotechnological production of aliphatic compounds.

High-resolution cryo-EM structures of respiratory complex I: Mechanism, assembly, and disease.

Respiratory complex I is a redox-driven proton pump, accounting for a large part of the electrochemical gradient that powers mitochondrial adenosine triphosphate synthesis. Complex I dysfunction is associated with severe human diseases. Assembly of the one-megadalton complex I in the inner mitochondrial membrane requires assembly factors and chaperones. We have determined the structure of complex I from the aerobic yeast Yarrowia lipolytica by electron cryo-microscopy at 3.2-Å resolution. A ubiquinone molecule was identified in the access path to the active site. The electron cryo-microscopy structure indicated an unusual lipid-protein arrangement at the junction of membrane and matrix arms that was confirmed by molecular simulations. The structure of a complex I mutant and an assembly intermediate provide detailed molecular insights into the cause of a hereditary complex I-linked disease and complex I assembly in the inner mitochondrial membrane.

Molybdate pumping into the molybdenum storage protein via an ATP-powered piercing mechanism.

The molybdenum storage protein (MoSto) deposits large amounts of molybdenum as polyoxomolybdate clusters in a heterohexameric (αß)3 cage-like protein complex under ATP consumption. Here, we suggest a unique mechanism for the ATP-powered molybdate pumping process based on X-ray crystallography, cryoelectron microscopy, hydrogen-deuterium exchange mass spectrometry, and mutational studies of MoSto from Azotobacter vinelandii. First, we show that molybdate, ATP, and Mg2+ consecutively bind into the open ATP-binding groove of the ß-subunit, which thereafter becomes tightly locked by fixing the previously disordered N-terminal arm of the α-subunit over the ß-ATP. Next, we propose a nucleophilic attack of molybdate onto the γ-phosphate of ß-ATP, analogous to the similar reaction of the structurally related UMP kinase. The formed instable phosphoric-molybdic anhydride becomes immediately hydrolyzed and, according to the current data, the released and accelerated molybdate is pressed through the cage wall, presumably by turning aside the Metß149 side chain. A structural comparison between MoSto and UMP kinase provides valuable insight into how an enzyme is converted into a molecular machine during evolution. The postulated direct conversion of chemical energy into kinetic energy via an activating molybdate kinase and an exothermic pyrophosphatase reaction to overcome a proteinous barrier represents a novelty in ATP-fueled biochemistry, because normally, ATP hydrolysis initiates large-scale conformational changes to drive a distant process.

A comparative study of single-particle cryo-EM with liquid-nitrogen and liquid-helium cooling.

Radiation damage is the most fundamental limitation for achieving high resolution in electron cryo-microscopy (cryo-EM) of biological samples. The effects of radiation damage are reduced by liquid-helium cooling, although the use of liquid helium is more challenging than that of liquid nitrogen. To date, the benefits of liquid-nitrogen and liquid-helium cooling for single-particle cryo-EM have not been compared quantitatively. With recent technical and computational advances in cryo-EM image recording and processing, such a comparison now seems timely. This study aims to evaluate the relative merits of liquid-helium cooling in present-day single-particle analysis, taking advantage of direct electron detectors. Two data sets for recombinant mouse heavy-chain apoferritin cooled with liquid-nitrogen or liquid-helium to 85 or 17 K were collected, processed and compared. No improvement in terms of resolution or Coulomb potential map quality was found for liquid-helium cooling. Interestingly, beam-induced motion was found to be significantly higher with liquid-helium cooling, especially within the most valuable first few frames of an exposure, thus counteracting any potential benefit of better cryoprotection that liquid-helium cooling may offer for single-particle cryo-EM.

Cryo-EM structure of respiratory complex I at work.

Mitochondrial complex I has a key role in cellular energy metabolism, generating a major portion of the proton motive force that drives aerobic ATP synthesis. The hydrophilic arm of the L-shaped ~1 MDa membrane protein complex transfers electrons from NADH to ubiquinone, providing the energy to drive proton pumping at distant sites in the membrane arm. The critical steps of energy conversion are associated with the redox chemistry of ubiquinone. We report the cryo-EM structure of complete mitochondrial complex I from the aerobic yeast Yarrowia lipolytica both in the deactive form and after capturing the enzyme during steady-state activity. The site of ubiquinone binding observed during turnover supports a two-state stabilization change mechanism for complex I.

Bayesian inference of rotor ring stoichiometry from electron microscopy images of archaeal ATP synthase.

The 'Bayesian inference of electron microscopy' (BioEM) framework makes it possible to determine the stoichiometry of protein complexes using 3D coarse-grained models and a relatively small number of cryo-electron microscopy images as input. We applied the method to determine the most probable rotor ring stoichiometry of the archaeal Na+ ATP synthase from Pyrococcus furiosus, a multisubunit complex able to produce ATP under extreme conditions. Archaeal ATP synthases consist of a catalytic A1 part and a membrane-embedded AO portion. The AO portion is composed of a rotor ring and the a-subunit. The rotor ring of P. furiosus ATP synthase is composed of 16-kDa c-subunits, each consisting of four helices forming a bundle, with only one Na+-binding site per bundle. This ring was proposed to be decameric from LILBID-MS analysis of the entire ATP synthase. By contrast, the BioEM posterior favors a c9 ring stoichiometry. With BioEM, we ranked coarse-grained models of the whole complex with different ring geometry, using 6400 unprocessed particle images of the A1AO complex collected in vitreous ice. BioEM makes it possible to probabilistically establish the domain stoichiometry using low-resolution information and comparably few particle images.

Structural basis for energy transduction by respiratory alternative complex III.

Electron transfer in respiratory chains generates the electrochemical potential that serves as energy source for the cell. Prokaryotes can use a wide range of electron donors and acceptors and may have alternative complexes performing the same catalytic reactions as the mitochondrial complexes. This is the case for the alternative complex III (ACIII), a quinol:cytochrome c/HiPIP oxidoreductase. In order to understand the catalytic mechanism of this respiratory enzyme, we determined the structure of ACIII from Rhodothermus marinus at 3.9 Å resolution by single-particle cryo-electron microscopy. ACIII presents a so-far unique structure, for which we establish the arrangement of the cofactors (four iron-sulfur clusters and six c-type hemes) and propose the location of the quinol-binding site and the presence of two putative proton pathways in the membrane. Altogether, this structure provides insights into a mechanism for energy transduction and introduces ACIII as a redox-driven proton pump.